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MASS SPECTROMETER INVESTIGATION THE ISOTOPES 
XENON AND KRYPTON RESULTING FROM THE 
FISSION U** THERMAL 


Abstract 


Mass spectrometer investigations have been made rare gas fission products 
extracted from uranium irradiated with thermal neutrons. The irradiated 
uranium rods were allowed stand for various periods time after irradiation 
permit the decay most fission product chains stable isotopes. Four 
stable isotopes xenon were found having mass numbers 131, 132, 134, and 136, 
and three stable isotopes krypton with mass numbers 83, 84, and 
the most abundant the latter group, probably formed directly fission. 
addition long lived krypton with mass was discovered which isomeric 
with 4.0 hr. reported previously. 

The relative abundances these isotopes which are related directly fission 
yields the corresponding mass chains have been determined with accuracy 
better. The mass numbers these fission chains can now identified 
with certainty comparing mass spectrometer abundance data with known 
yield values the active chain members. Finally, the half-life krypton 
was determined comparing its concentration that stable isotope over 
period time. 

Introduction 


the uranium fission process thermal neutron captured the 
This gives rise activated compound nucleus which 
turn splits into two fragments fission products addition several 
neutrons. The process can represented the following equation 


U235 + U236* a + v (n) 


The primary fission products and have more neutrons for their charge 
than the stable nuclei and undergo several steps ending 
stable nuclear structures. Some such fission chains are now 


Most the work with fission products has involved study the activity 
the separated active products. The difficulties involved getting clean- 
cut separations active products, determining mass numbers fission 
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product chains, and, finally, determining accurate fission yield data 
these methods, are well known. The availability Canada small samples 
rare gas mixtures (0.001 cc. gas N.T.P.) extracted from irradiated 
uranium for volume measurements and counting (2) experiments, made 
possible for the first time mass spectrometer study certain stable isotopes 
resulting from fission product chains. Because the large number stable 
isotopes xenon and krypton which block further fission 
product chains, the mass spectrometer experiments were expected give 
important information about relatively large number these chains. 

Some six rare gas mixtures extracted from irradiated uranium rods were 
investigated December 1945. These samples contained about cc. 
xenon and about cc. krypton. Altogether eight isotopes krypton 
and xenon were found, seven stable and one active. Since the rods were 
allowed stand for six nine months, only long lived active isotopes will 
present with the stable ones. Relative abundances these isotopes were 
obtained with accuracy better. Details the investigations are 
reported this paper. 

Experimental 
Rare Gas Samples 

The rare gases investigated were prepared Drs. Arrol, Chackett, and 
Epstein (2), who had apparatus set for the extraction and fractionation 
rare gas mixtures for volume determinations. Details the simple extraction 
are given their papers. uranium used was the form composite 
rod made disks varying weight. Two such rods were irradiated with 
neutrons pile. Rod was irradiated 1.5 times the neutron intensity 
Rod but for only two-thirds the exposure time. 

Preliminary mass investigations showed cu. mm. extracted 
gas contaminated with about 60% argon. The cu. mm. samples 
therefore contained less than 0.0004 cu. cm. krypton and xenon from 
uranium fission. The gas sample line (see Fig. was designed handle 
these small samples gas. consists capillary line cc. volume 
leading capillary leak which introduces gas samples directly into the 
ionization chamber the ion source system (see Fig. 2). This small volume 


connected sample tube manifold means glass break seals that 


can opened for pumping and for the introduction samples. addition, 
the line fitted with large reservoir mercury pump that used pump 
samples into the small volume behind the mass spectrometer leak. This 
pump also used mix the gas sample and regulate the pressure gas 
behind the leak, obtain suitable ion currents for identification 
isotopes and for abundance measurements. The small cc. volume can 
sealed off from the rest the line the second break seal manifold after the 
sample has been introduced into the large reservoir the mercury pump. 
After the investigation each sample, the residual gas can condensed 
charcoal liquid air trap and stored for further work. 
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TO VACUUM 


STANDARD SAMPLE MANIFOLD 


SAMPLE MANIFOLD 


—MANOMETER 


—ORY ICE 
TRAP 


PRESSURE 
REGULATOR, 


Mass spectrometer sample line for rare gas 


Mass Spectrometer 


The mass spectrometer used 180 degree Nier type described previously 
(6, ion source system and part analyser tube are shown 
Fig. Plate supporting plate and forms the top the ionization 
chamber. Plate contains the first collimating slit, mm. long and 0.35 mm. 
wide. Between Plates and located box through which electrons pass 
from filament through 0.2 mm. slit and ionization chamber collector 
box trap, Plate mounted the ionization chamber over the first 
collimating slit, made volts more positive than Plate order 
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draw the ions out the chamber. Plate contains the final collimating slit 
0.25 mm. width through which the ions enter the analyser tube. The 
tube, which included ion source assembly, analyser plate assembly, and 
180 degree 5-in. radius curvature analyser tube, completely enclosed 


SOURCE ASSEMBLY 


CAPILLARY 
Leak, 


FILAMENT 
SHIELD 


=, 


TO PuMPS 


Fic. Mass spectrometer tube. 


glass envelope. The glass envelope wrapped with heating element 
that the tube can thoroughly baked. vacuum pump leak connected 
directly the analyser tube ensure low pressure this region, and 
thereby reduce scattering the ion beams toa minimum. With this arrange- 
ment there exists pressure gradient between the source end and the analyser 
tube, which are connected means the slit Plate and between the 
ionization chamber where the gas samples enter and the source end the 
tube, since this chamber enclosed but for two small slits. 


d-c. amplifier with linear feedback used measure the ion currents. 
This unit, addition the electronic equipment for control electron 
emission, accelerating ions, and stabilization magnetic current, have been 
designed operate from standard a-c. power lines. Details these units 
and discussion their performances are described elsewhere (5). 
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Results and Discussion 

Mass Spectrograms 

Figs. and show typical mass spectrograms for normal and fission product 
xenon, and Figs. and for normal and fission product krypton, respectively. 
each case the ion current plotted against the ion accelerating voltage, 
which given terms the corresponding mass number. can seen 
from the figures, the mass peaks are completely resolved. check the 
abundances the xenon isotopes standard mixtures was made before the 
investigation fission product sample. each case the mass spectrometer 
ratios agreed with known values within the limits accuracy (1%). Fig. 
shows the existence four xenon isotopes the fission product sample, masses 
131, 132, 134, and 136, and Fig. shows the existence four isotopes 


POSITIVE 


Fic. mass spectrogram for normal xenon gas. 


krypton with masses 83, 84, 85, and 86. the identification these masses, 
all the usual precautions were taken and the possibility contamination was 
investigated thoroughly. the first sample investigated, xenon 129 was 
found present the extent part 2000 the total number 
xenon atoms. contamination with normal xenon was suspected 
because the high concentration naturally occurring xenon. Later 
work showed this the case. Contaminated samples were found have 
traces and proportion their abundances normal xenon. 
Other samples free contamination contained less than part 40,000 
these two isotopes, which the limit detection for the instrument with the 
size samples available. The largest samples investigated showed the 
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Fic. mass spectrogram for fission product xenon gas. 
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mass spectrogram for normal krypton gas. 


existence trace amounts krypton and 82. These traces amounted 
less than part 8000 and 5000, respectively, terms total krypton. 


Isotope Abundances and Fission Yields 


The samples investigated contained only from cc. 
xenon and krypton. was therefore necessary use very small volume 
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behind the mass spectrometer leak and operate relatively low pressures. 
Because this the ion currents decreased slowly with time. However, 
making time corrections, quite precise abundance data were obtained. The 


CURRENT 


POSITIVE 


MASS UNITS 


Fic. mass spectrogram for fission product krypton. 


ion currents for the different isotopes were compared. Tables and give 
typical abundance data obtained for the isotopes xenon and krypton 
mixture rare gas fission products (Rod 1). 


TABLE 


TYPICAL MASS SPECTROMETER ABUNDANCE DATA FOR ISOTOPES OF FISSION PRODUCT XENON* 


Mass unit 131 132 134 136 
Ion** Atom Ion Atom Ion Atom Ion Atom 
Set 
current per cent current per cent current per cent current per cent 

1 255 13.67 370.5 20.35 670 35.92 560 30.02 

2 241 13.76 358.9 20.52 633.1 36.18 516 29.49 

3 236.4 13.76 351.6 20.47 620.3 36.09 509 29.73 

4 1083.6 13.75 1602.4 20.33 2838.0 36.00 2355 29.87 

5 204.4 13.73 301.4 20.25 553.8 35.87 448 30.10 
Average 13.73 + 0.03 20.38 + 0.08 36.01 + 0.1 29.82+0.2 


Sample extracted 407 days after the irradiation period. 
Ion current given terms deflection centimetres. 


Tables III and give summary stable isotope abundances obtained 
with five samples rare gas mixtures. Samples were extracted from 
the centre disk single uranium rod referred Rod Sample was 
extracted from centre disk Rod which was irradiated times the 
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TABLE 
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TYPICAL MASS SPECTROMETER ABUNDANCE DATA FOR ISOTOPES OF FISSION PRODUCT KRYPTON* 


Mass unit 
Set Ion 
current 
1 215 
2 204.6 
3 186.5 
4 540.0 
Average 


Atom Ion 
per cent current 
“14.3 401 
14,25 383 
14.10 353 
13.82 1067 
14.12 + 0.15 


84 
Atom Ion 

per cent current 
26.61 106.2 
26.70 102.5 
26.75 93.0 
27.2 273.5 
26.81 + 0.2 


86 
Atom Ton Atom 
per cent current per cent 
7.06 781 51.90 
7.13 744 $1.75 
7.06 686 52.0 
7.0 2026 $1.8 
7.06 + 0.04 51.86+0.08 


Sample extracted for analysis 407 days after irradiation period. 
not stable isotope (half-life (see below). 


TABLE III 


SUMMARY OF ABUNDANCE DATA FOR ISOTOPES OF FISSION PRODUCT XENON 


Sample 


(Rod 


Average 
(Rod 2)* 


13.73 0.04 
13.38 0.04 


Mass unit 
132 134 
20.4 +0.1 36.0 +0. 
20.38 0.08 36.01 
20.09 0.04 35.76 +0. 


30.77 


Rod was irradiated 1.5 times the neutron intensity Rod but only for two-thirds 


the exposure time. 


TABLE 


SUMMARY OF ABUNDANCE DATA FOR STABLE ISOTOPES OF FISSION PRODUCT KRYPTON* 


Sample 


(Rod 


Average (1-3) 


(Rod 
(Rod 


Mass unit 
14.1 28.3 
15.8 29.2 
15.1 28.7 
15.23 0.15 28.90 0.2 
15.39 0.05 28.91 0.05 


57.6 
55.0 
56.2 


56.3 +1.2 


55.87 0.1 
55.70 0.05 


Active krypton not included since its concentration function time and would vary 


from sample sample. 


131 136 
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neutron density for two-thirds long irradiation time. The first three 
samples contained such small amounts krypton that was difficult 
the precision the abundance data. These three sets data have 
been averaged for comparison with results for Sample which 
ably larger and permitted estimation the precision. Because the 
higher concentration xenon the samples, about six times that for krypton, 
higher ion currents and therefore higher precision could obtained for the 


xenon isotope abundances. 


which appears along with and not included 
Table since radioactive and its concentration changes with time. 
each case the uranium rods were allowed stand for number months 
after irradiation permit the different fission product chains decay 
stable xenon and krypton isotopes. Samples extracted from Rod 
over period months, the first about nine months after irradiation, 
not show any variation the stable isotopic abundances within the limits 


our precision. 
TABLE 


COMPARISON OF XENON ISOTOPE RATIOS FOR GAS SAMPLES EXTRACTED FROM 


URANIUM RODS IRRADIATED UNDER SLIGHTLY 
DIFFERENT POWER CONDITIONS 


Mass ratio 


Rod No. 
132/131 134/131 136/131 
1.48 2.63 2.17 
1.50 2.67 2.30 
Difference, 1.3 
TABLE 


COMPARISON OF KRYPTON ISOTOPE RATIOS FOR GAS SAMPLES EXTRACTED 
FROM URANIUM RODS IRRADIATED UNDER SLIGHTLY 
DIFFERENT POWER CONDITIONS 


Mass ratio 
Rod No. 
84/83 86/83 
1.90 3.66 
1.88 3.62 
Difference, 1.0 


Tables and VI, the stable isotopic ratios for fission product xenon and 
krypton obtained for Rods and are compared. These results show that 
stable krypton isotope ratios for the two agree well within the limits our 


precision, but that the xenon isotope ratios differ for the two much 
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6%. would appear that even the small variations the 132/131 and 
134/131 ratios are significant, view the high internal precision the 
measurements. Certainly the large difference the 136/131 ratio signifi- 
cant and cannot explained the basis mass discriminations the 
mass spectrometer. this explanation were true, larger discrepancies would 
have shown for the krypton isotope ratios, which were examined the 
same time. Further, the variations indicated Table are not linear 
function the masses. 
Krypton 

Ion currents appeared corresponding mass all krypton fission 
product samples. Purification the krypton samples did not alter the results 
relative abundance measurements and there was evidence for contamina- 
tion the mass spectrometer tube the sample line that could give rise 
mass ion current. number samples have been investigated and 
now quite clear that product the fission process and radio- 
active with half-life about years. 


The discovery this isotope was little surprising view the existence 
four hour reported previously. The new isotope obviously 
isomer the four hour with long half-life. sample krypton 
extracted from uranium rod 407 days after its irradiation with neutrons 
formed pile was found contain 7.06 0.04% the active 
isotope (see Table concentration was found less older 
samples and greater younger samples, indicating definite activity. This 
activity was also investigated Arrol, Chackett, and Epstein, who pumped 
fission product krypton into gas counter (3). Recently the existence 
long lived krypton isotope was mentioned address its 
discovery being attributed Sugarman and Hoagland. Undoubtedly, 
this isotope, found independently counting experiments, the active 
positively identified mass spectrometer methods. 


The decrease the concentration with time relative one more 
the stable krypton isotopes was measured and its half-life calculated. The 
results are given Table these calculations the assumption made 
that variations the concentration relative the stable isotopes from 
sample sample are due wholly the decay with time, and not due 
variations the distribution fission product isotopes. This assumption 
seems justified for the first two calculations made, Rows and Table VII, 
since the samples investigated were extracted each case from the centre 
disks the same uranium rods, and the relative abundance measurements did 
not show any significant variations the distribution the stable isotopes 
krypton and xenon that result from fission. The third calculation made 
from data Row Table VII, perhaps less justified since the 144 day 
sample came from uranium rod different from that the 407 day one, and 
the relative abundance measurements for the two samples did show signifi- 
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TABLE VII 
HALF-LIFE FROM MASS SPECTROMETER ABUNDANCE DATA 


Time after irradiation, days half-life 
Ratio 
144 407 619 years 


Active 
This sample from uranium rod. No. all others from rod No. 


cant variations the concentration. However, the variations the 
relative concentrations the other stable isotopes xenon and krypton 
were small and perhaps not significant. 

The results the three half-life calculations agree fairly considering 
the small changes the concentration measured. The assumptions 
made above calculations become unnecessary when abundance measure- 
ments can made the same sample rare larger sample, which 
will sufficient for several investigations, and which can studied over 
many months, therefore being prepared. 

The fission yield the year can calculated terms from 
the data Table II. Considering the half-life years and the sample 
407 days old, the ratio 84/85 turns out 3.50. 


Ratio Xenon Krypton 

date has been difficult determine with high accuracy the ratio 
total xenon total krypton with the mass spectrometer, because the small 
size the rare gas samples available and the difficulties encountered cali- 
brating the instrument with standard rare gas samples that size. was 
found that, under the conditions necessary handle these small samples, 
there was considerable fractionation the sample the capillary leak. 
However, extensive calibration the instrument with small standard rare 
gas mixtures, preliminary value 6.2 has been obtained for the ratio 
xenon This compares with 5.0 obtained Arrol, Chackett, and 
Epstein (2), who fractionated rare gas fission product mixtures and measured 
the volumes the different components. knowledge this ratio together 
with the isotopic abundances makes possible tie fission yield data for 
mass chains ending xenon the heavy fission product group with those 
ending krypton the light group. 


Discussion 


The relative abundances the stable isotopes found fission material give 
the relative fission yields the different mass chains directly. view 
the direct character and accuracy the measurements, these fission yield 
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values should make possible draw the fission number curve 
with greater certainty. Actually fission yield data obtained counting 
experiments only good 10%. From knowledge the absolute 
fission yield the 131 mass chain obtained for the day fission yields for 
mass chains 132, 134, 136, and 83, 84, 86, follow directly from our data. 
This would give seven points the curve. Actually these points have been 
fitted the best available curvef, and, with the exception one point the 
values agree consistently. This one discrepancy decidedly significant and 
indicates some sort branching from one mass chain another. The 
problem being investigated further. 

The following series compiled from Seaborg’s table isotopes may lead 
stable xenon (8). 

Early branching 
ratio data (1) 


days 
min. 
min. days Stable 
min. hr. 2.4 hr. Stable 
min. min. min. Stable 
1.8 min. 


Seaborg his 1944 table isotopes refers the mass numbers for and 
131. Since only four isotopes xenon have been identified with the 
mass spectrometer, and definitely associated with the eight day 
the must associated with the series and above, 
respectively. The order seems fairly certain, since one might expect the 
iodines become less stable with increasing mass. Further fission yield 
data reported for the above series should agree with the relative abundances 
obtained for the corresponding xenon isotopes, and this way the mass 
numbers can assigned the various series with certainty. 

mentioned above, does not exist fission product xenon within 
detection limits part 40,000. the other hand, the fission yield the 
129 series determined from the activities and Fermi, and 
Grosse (1) This suggests the existence long lived which 
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has been suspected for some time since activity that could attributed 
has yet been found. 

The following fission product chains previously indicated (8) might give 
rise stable Kr. 


140 min. Stable 
min. Stable 
min. hr. 
(4) Br8?7? Kr87? 
sec. 


From the mass spectrometer results, can conclude that only four stable 
isotopes krypton exist fission material more than trace amounts, 
and The presence and was expected since 
the mass numbers for Series and above had been definitely established 
previously Se(d, and reactions, respectively (4, 9). 
However, previous mention had been made the possible production 
and was entirely unsuspected view the known four hour 
which would have completely decayed the month period that elapsed 
between the time irradiation and analysis. The assignment mass 
number Series considered certain Seaborg (8). The dis- 
covered with the mass spectrometer must therefore isomer the four 
hour with long half-life. This conclusion led the half-life deter- 
minations given Table VII. Further counting experiments Arrol, 
Chackett, and Epstein fission product krypton gas revealed active 
krypton with long period (3). 

There known evidence for which would parent 
Therefore either the parent has very short half-life and has escaped 
detection, formed directly fission, both. 


The mass spectrometer method particularly suited the study fission 
products. the only method that gives information about the stable 
isotopes that grow from active products and that may form directly fission. 
Mass numbers can positively identified and relative fission yield data for 
the stable isotopes give good deal information about chain relations and 
possible decay mechanisms. Further, active products with long half-lives 
can investigated and their periods determined. extracting fission 
product samples from irradiated uranium soon after irradiation, the method 
can extended the short period isotopes planned explore 
this method further with the investigation active rare gas samples. 
will considerable interest also follow changes the relative abundances 
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the different fission product isotopes function time after irradiation 
and the pile irradiation intensity. Indications are that 
significant variations occur the distribution fission products, depending 
the intensity irradiation (see Table V). Further, the accurate mass 
spectrometer determination the krypton isotope abundances furnishes 
very sensitive method for detecting changes fission product distribution 
effected changes the energy neutrons the irradiation process, since 
the krypton isotopes fall the steep part the fission curve. 
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MEASUREMENT THE DIFFUSION LENGTH THERMAL 
NEUTRONS GRAPHITE! 


Abstract 


The theory and method measuring the diffusion length thermal neutrons 
graphite are discussed detail. The graphite pile was rectangular paral- 
lelepiped, 185.8 cm. square and 153.6 cm. high. The source was 1.2 gm. 
radium mixed with beryllium. Part the pile was used slow down the 
neutrons thermal velocities, and the density distribution thermal neutrons 
was explored the remainder the pile with detectors dysprosium oxide. 
The length thermal neutrons this graphite average pile density 
1.55 gm. per cc. was found cm., with probable error cm. 


Introduction 


The diffusion length graphite required the design graphite— 
uranium piles, reflectors, and thermal columns. has also been used 
index the purity commercial graphites. the early work leading 
chain reacting piles the heterogeneous type was 
clear that only with some difficulty could the value the multiplication 
constant brought above unity and the neutron density the pile 
made diverge. the various steps taken increase one was increase 
the thermal utilization, one the factors making This was brought 
about improving the purities the uranium metal and the graphite 
moderator, thus minimizing the parasitic capture thermal neutrons (see 
Smyth report (4) 

This paper contains account the the diffusion length 
thermal neutrons commercial graphite produced Canada. The 
measurements were made the Montreal Laboratory during June and July, 
1943, and the final report was issued October 29, 1943. will made 
clear later that the accuracy the diffusion length depends great deal 
the size the graphite pile. had that time only nine tons which 
slow down the neutrons from mixed radium and beryllium source and 
measure their diffusion length while thermal, the accuracy was somewhat 
disappointing. 

Description Graphite Pile 


The graphite pile wcighed 9.1 short tons and was built unplaned bars 
The rectangular parallelepiped formed the layers bars each was 


Manuscript received August 22, 1946. 
Contribution from the Physics Branch, Laboratory, Atomic Energy Division, 
National Research Council Canada. Issued N.R.C. No. 1468. 
Physicist, National Research Council; now the Chalk River Laboratory. 
were aware the time this experiment that similar work had been done the United 
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185.8 cm. square and 153.6 cm. high. The average pile density was 1.55 
gm. per cc. while the densities planed bars taken random were found 
range from 1.58 1.65, with the average 1.62 The pile rested 
cadmium sheet, 0.3 mm. thick, covered with wrapping paper, which served 
define its lower boundary for thermal neutrons. 


PARAFFIN 


SOURCE 


Graphite pile. 


The neutron source, 1.2 gm. radium mixed with beryllium, was placed 
about cm. from the south face the pile. This face was covered with 
paraffin wax, about cm. thick, reduce the number neutrons escaping 
from and reaching other parts the pile scattering from the walls 
the room. 

Experimental Method 


The density distribution thermal neutrons the pile was measured with 
dysprosium detectors along the central axis, which the source 
was placed, and the vertical and horizontal lines that intercept this axis 
normally 13.5, 50.4, and 87.2 cm. from the north face. These lines may 
also described the defining lines three vertical planes, 
labelled first, second, and third planes order towards the north end the 
pile. The planed movable bars into which the detectors were placed are 
distinguished Fig. the shading their ends. The perpendicular 
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distances from the source the three vertical planes measurement and 
the north face were 81.3, 118.1, 155.0, and 168.5 cm., respectively. This 
disposition source and planes measurement was chosen give volume 
reasonable length (about 1.7 diffusion lengths) which the density distribu- 
tion thermal neutrons could explored, and the same time locate this 
volume sufficient distance from the source fast neutrons that the 
correction for the thermal sources within was small. The evaluation the 
diffusion length therefore did not require the precise value the slowing- 
down length the neutrons from the source. 

set detectors was used. consisted 280 mgm. dysprosium 
aluminium foil, were wedged holes the movable bars with half-cylindrical 
graphite plugs, and always faced south the pile. order avoid the 
influence one detector the activity acquired another, adjacent ones 
given time the pile were never less than cm. apart. 

The relative sensitivities the detectors were determined exposing 
them turn neutrons under constant geometrical conditions the pile 
and comparing their with that from radium preparation. 
Such standards were used systematically throughout the measurements 
minimize errors due variations the behaviour the 
counters and associated amplifiers and scales-of-32. 

After the activity had approximately reached saturation point the 
pile, eight minutes from the instant removal was allowed elapse for short- 
period activity vanish, and the dysprosium activity was compared with 
radium standard bracketing for about two hours. The density 
thermal neutrons was (arbitrarily) expressed the counting rate 
per minute), corrected saturation, the instant removal 
from the pile, taking 140 min. for the half-period dysprosium. Corrections 
were also applied for any activity remaining from the previous activation, 
and for the long period activity (perhaps contamination) the detector. 
The latter was from two four per minute from single detector. 
The half-period dysprosium was checked about comparing decay- 
curves with those actinium active deposit, the which (36.1 
min.) accurately known. 

Corrections had applied the measured distribution neutrons 
allow for two extraneous effects, which would have been greatly reduced 
had had sufficient cadmium sheet and covered the pile. (i) Some neutrons 
entered the pile from external sources, used other experiments the same 
room. (ii) Some neutrons from the source the pile were reflected from the 
walls the room back into the pile. The contribution these neutrons 
the measurements was investigated with the dysprosium detectors measur- 
ing the densities few points the pile the absence external sources, 
and many points the pile the presence external sources and 
with our source placed paraffin wax box the floor beyond the south 
face. arrangement our source was retained the room reproduce, 
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least part, the effect scattered neutrons from the walls the room. 
The final correction, representing the effects all neutrons from outside the 
pile, course varied from point point the pile but never exceeded nine 
per minute. When these intensities were subtracted from the 
original intensities, the extrapolated distributions vanished the theoretical 
bounding planes for such pile vacuo. 


shall now examine our fundamental assumption that the density 
thermal neutrons proportional the saturation activity our detectors 
facing south the pile. order have reasonably high counting rates 
were obliged choose rather thick detectors having high capture cross- 
section per atom. When stronger source became available, the 280 mgm. 
detectors, orientated several directions, were compared with mgm. 
detectors dysprosium oxide number points the north-south axis 
the pile. The 280 mgm. detectors acquired about more activity when 
facing south and less activity when facing north than the average for 
different points between the first and third planes. The 
mgm. detectors showed measurable orientation effect. The ratio the 
averaged activities acquired thick and thin detectors was constant 
4.6 within experimental error. The 280 mgm. detectors, used this 
experiment, therefore give sufficiently accurate measurements the density 
thermal neutrons. small departure from the law capture 
neutrons dysprosium, suspected some experimenters, can scarcely 
importance here. 


Elementary Diffusion Theory 


The diffusion equation for thermal neutrons may simply derived from 
the idea continuity the steady state: 


production capture outflow, 


the number thermal neutrons produced slowing down per cm? per sec. 
the density number thermal neutrons per cm? 

the mean lifetime with respect capture. 

the current net number per cm? per sec. 


analogy with the kinetic theory gas diffusion introduce the 
diffusion coefficient and the assumption grad the equation 
continuity becomes: 


The product Dr, having the dimensions the square length, written 
where defined the diffusion length. Equation (2) becomes: 
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assume that there are fast neutrons the region beyond the first 
plane, which 81.3 cm. from the source, Equation (3) becomes: 


(4) 


order proceed find function the co-ordinates, the origin 
chosen the lower north-east corner the pile, with the positive direction 
the X-axis westward, the Y-axis upward, and the Z-axis southward into 
the pile. The solution the diffusion Equation (4) subject the boundary 
conditions: 


where 189.4 cm., 157.2 cm. are the two dimensions the pile 
corrected for the breakdown the diffusion theory near the boundaries. 
This correction amounts shift the theoretical boundary vacuum 
distance 0.71 times the transport mean free path. This theoretical result 
was brought our attention Placzek and Seidel’s work problem 
transport theory (3) (see also Hopf (1)). The remaining boundary 
condition given the densities measured the first plane, which may 
analysed into Fourier components: 


The distribution assumed symmetrical about the axis 
the pile, which was confirmed experimentally. The required solution 


Equation (4) is: 


sinh 

lm a b inh d 

Lit 

l,m a b 


The are called the relaxation lengths the pile for the various com- 
ponents. The analysis the measurements aims isolating and evaluating 
and turn the diffusion length through Equation (7). 


Experimental Results 


Fig. shows the densities thermal neutrons plotted against distances 
measured vertically and horizontally from the north-south axis the pile. 
The density each point was measured twice, and, taking corresponding 
points opposite sides the axis together, each density plotted Fig. 


| 
| 
| 


CANADIAN JOURNAL RESEARCH. VOL. 25, SEC. 


the average four measurements. The ordinates are the counting rates 
per minute. For points the first and second planes, where 
the counting rates are reasonably high, separate determinations agreed 
about 3%. 


2600 


DENSITY THERMAL NEUTRONS 


DISTANCE FROM CENTRAL AXIS PILE 


Fic. Measured densities thermal neutrons along horizontal lines three vertical 
nes are shown circles, and those along vertical lines the same planes are shown 
filled circles. The curves shown are those for best fit according the method least squares. 


The amplitudes were found from the amplitudes the trigonometric 
functions required fit the experimental data along each axis each plane. 
The method least squares was used. Neglecting and higher Fourier 
terms, the following coefficients were obtained. 


First plane 89.0 cm.) 


Aun = 2290, An => 166.5, Ais => 68.1 (An = 26.1, Ais os 28.4). 
Second plane 52.2 cm.) 


Third plane 15.3 cm.) 
The fit these values shown the curves drawn Fig. 
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From Equation (5) 


sinh 89.0 
2290 
7 
Lu 
Hence 
Similarly, 
135.6 15.3 
Hence 


The latter values are less reliable than the former the difficulty 
measuring the small intensities the third plane. correction forthermal 
sources near the first plane must applied, and this matter will dealt 
with the next section. 


The densities measured along the central north-south axis the pile are 
shown Fig. point represents the average least two measure- 
ments. few points were remeasured later with stronger source, and these 
results are distinguished filled circles. These activities, averaged for 
different orientations the detectors, were normalized fit the earlier 
measurements the first point. The curves were calculated for 48,52, 
and cm. and slowing-down length 20.8 cm. 


Corrections for Sources Thermal Neutrons 


Consideration the effect thermal sources the derived diffusion length 
will based the slowing-down theory. convenient divide the 
discussion into three parts. 


(1) Effect the Paraffin Wax the South Face 


The wax did not increase the number epithermal neutrons beyond the 
first plane measurement any extent, for the total distance from the 
radium source the wax and back the first plane 117 cm. more than 
five slowing-down lengths. The wax may have increased the number 
thermal neutrons the region measurement, compared with the case where 
the south face boundary air. This point was tested experimentally, 
and was found that the removal the wax caused maximum reduction 
density 10% the first plane. 


The experimental arrangement with the wax reflector was closer the 
case square pile, bounded air all sides, than the case semi- 
infinite pile the z-direction, for calculation shows that the density thermal 
neutrons the first plane would reduced about 40% going from the 
semi-infinite the square pile. The effect the wax was therefore reduce 
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the ratio fast slow neutrons the region measurement 10%; 
and hence slightly smaller correction for thermal sources required than 
the case the finite pile. 


Z 


THERMAL NEUTRONS 


CM. 


DISTANCE ALONG AXIS PILE 


Fic. densities thermal neutrons along the central axis the pile are 
shown unfilled circles. using stronger source and after normalization with 
the measurements the first plane are shown filled circles. The lowest curves are 
those calculated for 48, 52, and cm. and slowing-down length 20.8 cm. 
addition, the curve shown for (no capture), and the straight lane for pile infinite 


being used the text for two the pile dimensions. 


(2) Correction for Thermal Sources the Case the Finite Pile the Basis 
Elementary Slowing-down Theory 
The equation for the number neutrons per unit volume passing through 
t 


0 

(change with age influx capture source), where Dirac’s delta- 


age which neutrons become thermal and liable captured. the 
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slowing-down length thermal velocities. the strength the source 
q 


pile the z-direction. 

take, first, the case pile infinite the z-direction and later use 

This solution, representing the distribution sources thermal neutrons, 
used boundary condition solving Equation (8) for thermal neutrons 


Since x(@) the number neutrons per unit volume passing through 
age time dt, the number neutrons per unit volume age interval 


given 


and the total density thermal neutrons given where 
the diffusion coefficient thermal neutrons. 


1 
where 


1 ao 


~(@,-2) 


equivalent treatment has been given Wallace and LeCaine (5). 


The boundary conditions the are now introduced super- 
posing the solutions for image sources about the two boundaries, viz., for 
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For the region cm., between the first and third planes, need 
only consider the primary source and three images. The first Fourier com- 
ponent the density for the four sources the only term for which the 
correction for thermal sources may significant fraction the total 
density. Neglecting contributions less than 0.1%, reduces the 
following: 


-% 2-2 


~(2;-2) 


estimate the slowing-down length may made from measure- 
ments taken Dr. Elliott. used indium resonance detectors and 
the same source this pile. sufficient for the present purpose state 
that single Gaussian curve fits the measured distribution near the source 
the slowing-down length indium resonance (1.44 ev. (2) taken between 
single Gaussian curve fits the measured points farther from 
the source the slowing-down length taken cm. Allowing for 
slowing-down from the indium resonance region the region capture 
dysprosium, was taken between and cm. 

Equation (9) for differs from Equation (5), where the variation 

find that these two terms amount —2.2% for 20.8 
cm., —1.5% for cm., and —0.7% for 18.5 cm. the first 
plane, and become negligibly small for distances from the source greater than 
100 cm., i.e., the region the second plane. 


The curves Fig. were computed using 20.8 cm. Equation (9) 
for pu, and the higher components from Equation (5). 


(3) Modification the Slowing-down Theory 


Part (2) the contribution thermal sources the density was calculated 
the basis elementary slowing-down theory. Using better approxima- 
tion the transport theory, for which are indebted Dr. Marshak, 
found that the additional correction the density distances 
100 cm. from the source never amounts more than 25% the correction 
found Part (2). 

The final correction the density for thermal sources approximately 
the first plane, and increased this amount, the value 
calculated from the ratio reduced cm. 


where 


with given sinh the last two terms the bracket. Using 
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Fig. the calculated curves are normalized the density the first 
plane. The best value fit the points the central axis between 
and cm. from the first plane (i.e., the region where the sensitivity 
greatest) cm., with probable error about cm. 

For this graphite average pile density 1.55 gm. per cc. the diffusion length 


Limitation Accuracy 


The density measurements were limited accuracy the rather low 
counting rates, instability the counter sensitivity, and the unsatisfactory 
reproducibility geometry the counter. 

The size the graphite pile was insufficient for good accuracy, which may 
seen from Fig. The straight line shown this figure represents the 
variation density along the axis pile infinite the and y-directions 
(no leakage neutrons through its sides) and (no capture). The 
next curve refers the experimental pile regards dimensions (same leakage) 
but having (no capture). clear that for pile with this cross- 
section the leakage neutrons more important than the capture within 
the graphite. The same point may shown algebraically. The relaxation 
length the first component the harmonic analysis given 


The second leakage term 1.8 times the first capture term. 
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THE DIFFUSION LENGTH THERMAL NEUTRONS HEAVY 
WATER CONTAINING LITHIUM CARBONATE! 


Abstract 


The density distribution thermal neutrons was measured with small boron 
trifluoride chamber cylindrical tank containing 113 litres heavy water 
which lithium carbonate was dissolved. The diffusion length was found 
22.7 cm. this solution containing 7.70 atoms lithium per molecule 
heavy water (99.4 atom After corrections were applied for the 
capture neutrons the heavy water and light hydrogen, the capture cross- 
section lithium was found cm.? per atom for neutrons 
standard velocity 2200 per sec. from the measured diffusion length and known 
transport mean free path. 


Introduction 


The diffusion length thermal neutrons solution uranium com- 
pound heavy water not easy measure directly owing multiplication. 
the ratio the capture cross-section uranium that boron lithium 
known, measurement the diffusion length solution boron 
lithium compound leads the diffusion length uranium solution. 


While the cross-section boron perhaps better known than that 
lithium, there strong reason against its use and favour the use 
lithium. Not only must the compound dissolved heavy water 
free ordinary water and impurities that capture thermal neutrons strongly, 
but must possible remove completely from the heavy water 
distillation. Mr. Hébert and Dr. Maddock (7) have kindly investi- 
gated these points for us, and found that the distillates from borax solutions 
contain boron, but that the lithium carbonate available this laboratory 
meets all the requirements. 


The density distribution thermal neutrons was explored with small 
boron trifluoride chamber cylindrical tank containing 113 litres heavy 
water which was dissolved sufficient lithium carbonate reduce the 
diffusion length suitable value for measurement. this solution contain- 
ing 7.70 atoms lithium per molecule water (99.4 atom D), 
the diffusion length was found 22.7 cm. This experiment was done 
the Montreal Laboratory during the winter 1944 and described 
report bearing the date May 16, 1944. 
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Theory the Experiment 
The diffusion equation the absence thermal sources 


where the density thermal function the co-ordinates, 
and the diffusion length (8). Equation (1) written cylindrical co- 

The origin taken the centre horizontal cross-section near the bottom 
the tank, and the positive direction the z-axis upwards along the axis 
the tank. The term involving the angular co-ordinate omitted from 
Equation Any angular asymmetry found experimentally can allowed 
for taking the angular average the measured density for given radial 
distance and given height 

The solution Equation (2) subject four boundary conditions. 

(1) The density finite the axis the cylinder. 

(2) The density vanishes the surface vertical coaxial cylinder 
radius where (evaluated later) the effective radius the tank; 


(3) The density vanishes horizontal plane short distance above the 
surface the liquid; 


(4) The density distribution the horizontal plane through the origin 
given measurement. 
The desired solution Equation (2) 


sinh 
= 1,2,3,... sinh 
where 
(4) 


The are the relaxation lengths the various components. 
The can found from the roots the equation: 


2.4048 5.5201 


Boundary condition (4) may written 
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This horizontal plane containing the origin must taken few scattering 
mean free paths above the bottom the tank avoid effects scattering 
the materials below the heavy water. 

The density distribution was measured three horizontal planes spaced 
12-cm. intervals, and along the central axis. The lowest first plane was 
approximately 4.2 cm. above the bottom the tank. The distributions 
were analysed transversely into Fourier—Bessel series, and the first component 
was thereby isolated. Its diminution the vertical direction permitted 
and found according Equation (3). The diffusion length was 
derived from Equation (4). 


Experimental Details 


The source thermal neutrons was pile 186 cm. 186 cm. 
horizontal section and 274 cm. high. the vertical axis the pile and 
100 cm. from the top was placed 2.2 gm. radium mixed with beryllium. 
Since the distance from this source fast neutrons the cylindrical tank 
top the pile was five slowing-down lengths, neutrons entering the tank 
were almost entirely thermal 99% number). 

The graphite pile was covered its sides, top, and bottom with cadmium 
sheet, 0.4 mm. thick. wooden platform, 3.5 cm. thick, rested the 
cadmium top the pile. Circular holes, cm. diameter, the 
cadmium top and the platform, were centred roughly the axis the 
pile. Acircular aluminium tray, cm. diameter, was placed symmetrically 
over these holes catch any leakage heavy water case accidents, after 
the space between the graphite blocks and the upper surface the wooden 
platform had been packed with petroleum coke and levelled. The experi- 
mental tank (Fig. was centred the tray. This tank was made 
aluminium wall thickness 3.3 mm. and 45.4cm. high. Its internal diameter 
was measured several directions and various heights with suitable 
inside caliper made for the purpose. The individual measurements were 
between 61.16 and 60.94 cm., and the average diameter was 61.01 cm. 

The cover the tank was designed with the assistance Dr. Bauer. 
glass disk, mm. thick and 69.7 cm. diameter, rested the upper 
flange the tank, and the joint was sealed with Apiezon grease radial 
slot, cm. long and 2.6 cm. wide, along which the boron chamber could 
moved, was cut the disk. glass slab, cm., with identical 
slot, was cemented the disk that the two slots coincided. second 
slab, cm., with conical hole for the gland supporting the boron 
chamber, rested the upper greased surface the first slab. tall box 
housing the head amplifier and gland for the boron chamber was clamped 
the movable sides the tank were tightly wrapped with cadmium 
sheets, and similar sheets were placed top the glass disk and slabs. 

number excellent boron chambers were made for Mr. Veall 
(11); one these shown Fig. chambers were made lead glass, 
mm. internal diameter and 0.45 cc. volume, and were filled with boron 
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trifluoride pressure about atm. The guard-rings were grounded 
the solution. The first chamber was sealed the lower end single glass 
tube housing the lead the central anode. The high tension the 
platinum cathode was rubber-insulated wire the lower end the chamber 
and the joint was made water-tight with rubber Bostik. After some time the 


Fic. Experimental tank. The boron chamber and head amplifier box are shown. 


insulation failed, and this chamber was replaced second chamber sealed 
two parallel glass tubes, one tube housing the lead the anode, and the 
other the high tension wire the cathode (Fig. 2). The latter tube was 
filled with some the solution reduce the disturbance the neutron 
distribution caused the channel. The glass tubes the chamber were 
clamped gland, and the height was adjustable. The gland for the second 
boron chamber, designed Mr. Lawrence and Mr. Veall, 

The counting equipment was supplied Mr. Freundlich. 
consisted head amplifier (5), mounted closely possible the boron 
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chamber, linear amplifier four stages, discriminator reject pulses 
below chosen size, scale-of-32, and oscilloscope. stabilized high 


CM. 


Fic. support boron chamber. Tightening the nut pushes the brass piece 
the rubber gasket thus making tight joint. 


Att. 
Fic. Small boron trifluoride chamber. 
D 5 
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voltage (usually about 800 v.) was applied the boron chambers. 
electrical leads were shielded, and all units were brought common ground. 

The solution lithium carbonate was made stirring 179.73 the 
carbonate 125.2 0.1 kgm. heavy water for period three hours. 
portion (125.1 kgm.) the solution was placed the experimental tank. 
After the experiment had been completed, Mr. Hébert analysed two samples 
the solution for and found 1.566 and 1.575 gm. lithium carbonate 
per litre, which may compared with 1.585 gm. per litre from the original 
weights the components. 

The cylindrical co-ordinates the approximate centre sensitivity the 
boron chamber the tank were specified follows (Fig. The angular 
co-ordinate was read the axial slot the circular scale painted around 
the tank. The radial distance was found from the displacement the 
movable glass plate supporting the boron chamber, the glass disk being 
centred through measurements between its edge and the upper flange the 
tank. The height the boron chamber above the first plane was measured 
from mark the gland mark the glass tube. 


Experimental Results 


any point the counting was carried out two parts, (i) with the 
boron chamber and gland given position, and (ii) with the boron chamber 
and gland rotated through 180° from the first position. average the 
two readings counts per minute was taken measure the density 
thermal neutrons the point. This method eliminated possible error due 
the fact that the chamber was not situated accurately below the centre 
the gland. Sufficient particles were counted each point for the statistical 
error less than 

Most our measurements the three planes are summarized Table 
the first plane, readings were taken along eight radii spaced 45° intervals. 
That is, the densities values 10, 15, 17.5 etc. were 
compared with the density for eight different angles readings 
were normalized count the centre 700 particles per minute, which 
was approximately the observed counting rate. Similar measurements were 
made the second plane, which was cm. above the first plane. this 
case readings were normalized count the centre 270 particles per 
minute (see Table slight angular asymmetry appeared both planes, 
and was traced variations depth due the fact that the bottom the 
tank was not perfectly flat and also slight error levelling. The 
asymmetry was removed averaging the densities for given and over 
the angle the third plane, cm. above the second, was sufficient 
measure along four radii. These readings were normalized count 
the centre particles per minute (see Table II). The average densities 
are plotted against the radial distances Fig. 
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TABLE 
EXPERIMENTAL RESULTS 
DENSITY VARIOUS ANGLES 
45° 90° 135° 180° 225° 270° *Av. 


700 700 700 700 700 700 700 
637 640 649 652 636 641 641 642 
567 568 568 585 578 571 574 560 571 
458 470 481 472 460 460 463 465 
392 393 406 405 390 394 401 
320 331 330 328 322 307 308 320 
Second plane 
270 270 270 270 270 270 270 270 270 
235 240 241 238 244 242 236 235 239 
197 208 207 207 206 202 207 205 
184 181 181 181 181 178 181 
153 161 160 158 157 151 157 
Third plane 
85.5 89.5 87.5 86.7 87.3 
80.5 78.1 79.0 79.6 79.3 
67.0 68.6 65.9 65.4 66.7 
17.5 59.5 58.5 60.5 59.5 
48.5 50.6 49.0 49.4 49.4 
22.5 39.5 39.8 39.6 
30.3 29.4 29.8 29.8 
27.5 18.0 18.0 


the last line 27.5 cm.), where there single measurement 270°, the 
average over all angular positions takes into account the angular asymmetry this position judged 
from the measurements neighbouring points and the contour map the bottom the tank. 


The density distribution along the central axis was also measured with the 
two boron chambers, and the averages points spaced 2-cm. intervals are 
plotted Fig. densities the centres the second and third planes 
were each compared many times with the density the centre the first 
plane. These results including the averages the residuals are given 


First plane 
i 
t 
i 


HEREWARD AL.: DIFFUSION LENGTH THERMAL NEUTRONS 


TABLE 
COMPARISON THE DENSITIES THERMAL NEUTRONS THE CENTRES THREE PLANES 


Boron chamber end Boron chamber end 
vertical tubes L-shaped tubes 
Plane No. 
Density comparisons Density comparisons 


700 


600 


FIRST PLANE 


DENSITY THERMAL NEUTRONS 
a> 


100 


RADIAL DISTANCE (CM) 
Fic. radial distribution the density thermal neutrons. 


All measurements taken with the two different boron chambers agreed 

closely and were independent the bias-setting the discriminator and 
the potential difference across the boron chambers. order ascertain 
whether there was any disturbance the neutron density caused the 
hollow glass tube supporting the boron chamber, the tubes the second 
boron chamber were bent into L-shape that the centre sensitivity 
the chamber was displaced 10.2 cm. horizontally from the vertical supporting 
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tubes, and the intensities the centres the three planes were again com- 
pared. From these results Column with those Column 
Table clear that the channel effect small. 


800 


700 


a 


500 


400 


DENSITY THERMAL NEUTRONS 


200 


DISTANCE FROM LOWEST PLANE 


Fic. The axial distribution the density thermal neutrons, unfilled circles represent- 
ing measurements with the boron chambers and filled circles those with dysprosium detectors. 
The position the liquid surface indicated the line segment 34.6 cm. 


few measurements along the central axis were made with dysprosium 
detectors. Each detector consisted 280 mgm. dysprosium oxide 
aluminium tray, 2.3 Several these detectors were spaced 
thin strip aluminium, which was slid into aluminium sheath oblong 
cross-section and wall thickness mils. This sheath was supported its 
top cone that fitted into the hole the movable glass slab. Our counting 
procedure for activated detectors has been described previous paper (8). 
The measurements with the dysprosium detectors are fitted the measure- 
ments with the boron chambers Fig. where the agreement seen 


satisfactory. 
Analysis the Experimental Results 


Before the transverse measurements could analysed into Fourier—Bessel 
series the parameters had estimated. This step required the value 
the effective radius our first analysis, was taken equal 
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the internal radius the tank, plus 0.05 cm. for the aluminium wall actual 
thickness 0.33 cm., plus 1.65 cm. This last length was the intercept the 
curve, showing the density thermal neutrons sample the same heavy 
water function the distance from flat cadmium boundary, obtained 
experiment Auger, Munn, and Pontecorvo (1). 


30.50 0.05 1.65 32.2 cm. 


The density distribution each plane may written 


(6) 
where 
sinh 
sinh 


The experimental values and were substituted Equation (6) 
form set experimental equations for the amplitudes 
method least squares was applied these equations. The amplitudes, 
derived from the normal equations, are shown Table III. The curves 
shown Fig. were drawn from Equation (6) with these amplitudes. 


TABLE III 


FIRST ANALYSIS 


(s) 814 287 
As (z) = 8 


The reduction the amplitude the first component the 
analysis with increasing used find the relaxation length. When the 
three planes are equally spaced, two useful formulae can derived from 
the first component the planes having vertical co-ordinates and 
2c, respectively. easy show that 


’ q 
§ q 


CANADIAN JOURNAL RESEARCH. VOL. 25, SEC. 


Similar expressions may written for other components. From Equations 
(7) and (8) and the amplitudes (before rounding off their values for Table III) 


11.63 cm. 
37.6 cm. 
From Equation (4) 23.5 cm. 


the corresponding equations are applied the second and third com- 
ponents, find that imaginary. Two possible explanations for this 
behaviour may suggested. (i) Since the relaxation length for the second 
component only 5.7 cm., compared with 2.3 cm. for the transport mean 
free path, the diffusion theory may greatly error for components higher 
than the first. (ii) Since 1/L? only (but 33% the deter- 
mination from the second (and higher components) too dependent 
and hence the effective radius. there some doubt regarding the 
contribution the aluminium wall the tank the effective radius, value 
32.5 cm. for the radius admissible. Using 32.5 cm. and 


the amplitudes (z) Table were obtained the method least 
squares. From Equations (7), (8), and (4) 


TABLE 


SECOND ANALYSIS 


First plane Second plane Third plane 

As (z) = 4 = tial 


and the amplitudes (before rounding off their values for Table IV) 


11.62 cm. 
cm. 
b, = §.78 cm. 
30cm. 


would possible derive perfectly consistent set values from 
the first and second components choosing suitable value very close 
32.5 cm. 

comparison the analyses for radii 32.2 and 32.5 cm. with the experi- 
mental results Table shows that impossible choose between them 
the basis best fit. However, change the effective radius 
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TABLE 


EXPERIMENTAL RESULTS AND ANALYSES 


700 700 700 

686 684 685 
642 643 643 
571 571 571 
17.5 525 523 523 
465 466 467 
22.5 401 401 400 
320 321 319 
222 222 222 

Second plane 

270 270 269 

261 262 262 
239 239 240 
205 204 204 
181 182 181 
157 156 155 
22.5 126 127 127 
96.0 95.4 95.8 
62.5 62.1 63.7 

Third plane 

91.0 91.0 91.0 
87.3 88.3 88.2 
79.3 80.2 79.9 
66.7 67.0 66.6 
59.5 58.7 58.4 
49.4 49.6 49.4 
22.5 39.6 39.7 39.7 
29.8 29.4 29.6 
18.0 18.9 19.5 


causes only change derived from the first component. have 
information the breakdown diffusion theory for components 
higher than the first, which was mentioned possible explanation the 
internal inconsistency found the first analysis. Accordingly, slightly 
prefer the second analysis with its better internal consistency, and shall write 


for the final 
or e final resu L = 22.7 cm. 


The axial measurements support, but not throw new light on, the value 
the diffusion length just derived. Fig. the measurements taken 
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2-cm. intervals are shown good agreement with the curve plotted from 
(3) with the following constants: 


While the liquid surface shown Fig. 34.6 cm. above the first plane, 
this value can only approximate, the centre sensitivity the boron 
chamber was not accurately located. (The chamber counts over length 
between 2.0 and 2.6 cm.) Owing this reason and the scattering 
neutrons from the glass cover not possible evaluate from the position 
the liquid surface and the value the cut-off distance the absence 
From inspection Fig. must lie between 36.0 and 
38.5 cm. However, the diffusion length determined mainly from the 
measurements deep within the and therefore not very much 
cependent 

The effect errors the observations the derived diffusion length can 
seen from the following considerations. The density distribution 
thermal neutrons obeys the differential Equation (2), which may written: 


With reference Fig. found trial that the curve altered 


Fig. gives poor fit with the points. The two right-hand terms 
(9) are roughly the ratio better, except near the wall 
the tank where less weight attached the observations. From the 
estimates the resulting error cannot exceed 


Other errors, such that the concentration lithium, and those due 
the channel effect and shift the centre sensitivity the boron chamber 
raised the solution, are comparatively negligible (Table 
calculation found that the shift the centre sensitivity the chamber 
causes change the counting rate less than 0.1%, which beyond 
detection. 


4-3 


Discussion 


the basis diffusion theory possible calculate any one the 
three quantities diffusion length transport mean free path and capture 
mean free path provided that the others are given, for 


(10) 


x 
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Auger, Munn, and Pontecorvo (1) have found 2.31 cm. sample 
the same heavy water used here for the measurement the diffusion length. 
Since the amount lithium carbonate dissolved the water too small 
affect its scattering properties, 2.31 cm. the solution. Substituting 
22.7 cm. and 2.31 cm. Equation (10), 


the number atoms (or molecules) per cubic centimetre the constituent, 
and the capture cross-section per atom (or per molecule), the sum being 
taken for the constituents and the solution. order 
derive the capture cross-section lithium, shall allow for the capture 
neutrons the heavy and light water,* and neglect the capture the 
carbon and oxygen the carbonate and the oxygen the light water present. 


based mass spectroscopic analyses carried out Dr. Thode’s labora- 
tory. For lithium carbonate, 1.57 gm. per litre solution, 


2.58 atoms lithium per cc. 


According argument given later, the capture cross-sections 
used here should those for neutrons 2500 per sec. Booker, 
Cavanagh, Hereward, Niemi, and Sargent (2) have found from measurement 
the diffusion length heavy water that the capture cross-section 
0.81 per molecule for neutrons velocity 2500 per sec. 

The capture cross-section hydrogen will found from the absolute cross- 
section boron and the ratio the cross-sections boron and hydrogen. 
Dr. Fermi (4) informed that the capture cross-section boron 705 
per atom for neutrons velocity 2200 per sec. (Recently Marshall 
(9) has published his value 710 per atom, measured some time 
ago the Argonne Laboratory.) The best measurement known 
2270, obtained Whitehouse and Graham (12) the Montreal 
Laboratory. From these results 0.548 per molecule 
for neutrons velocity 2500 per sec. 

for the lithium. 
The capture cross-section lithium follows: 


For cross-section calculations immaterial whether consider 1.2% HDO its equiva- 


lent 0.6% 
tWe not make the classic distinction between velocity and speed. 


| 
2 
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will seen from the considerations that follow that the cross-section 
found from the diffusion length and transport mean free path corresponds 
approximately the average velocity the neutron spectrum, which 
about 2500 per sec. room temperature. assume that the transport 
mean free path heavy water independent the neutron velocity, and 
therefore the diffusion coefficient 


D = 3h, 


where the average speed the Maxwellian distribution neutrons 
room temperature. The diffusion length may written terms the 
diffusion coefficient and the mean lifetime 


Since the capture cross-section lithium varies inversely the velocity, 
according the measurements Havens and Rainwater (6, 10), 
independent the velocity. Therefore 


where the capture mean free path for the average velocity the 
neutron spectrum, and hence found above the capture cross-section 


corresponding that velocity. Referring standard velocity 2200 
per sec., have 


view the assumptions and probable errors the combined physical 
data just used, conclude that this value agreement with 
Fenning, Graham, and Seligman (3) from their measured ratio and 
the Argonne Laboratory value 705 other hand, 
Havens and Rainwater (6, 10) have found higher values: 740 
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THE CALCULATION LINE STRENGTHS FROM 
LABORATORY DATA! 


PETRIE? 


Abstract 


King has made extensive intensity estimates spectral lines furnace 
and arc spectra, and these data provide useful material for determining line 
strengths. The estimates have been made for many elements astrophysical 
interest, and are reduced line strengths for Fel, Til, and Till. 


Introduction 


well distinguish between two quantities, the strength, and the 
intensity, spectral line. The strength line defined the sum 
the absolute squares the matrix elements electric moment joining the 
sets states question. the other hand, the intensity line propor- 
tional the number atoms any one the states its initial level, the 
fourth power its frequency, and its strength. 

The line strength the transition from initial final level defined 
the quantum numbers yLSJ and respectively, 


where the orbital quantum number, 
the spin quantum number, 
the inner quantum number, 

the magnetic quantum number, 

Expression (1) enables one calculate the relative strengths the lines 
within multiplet when coupling holds. Such line strengths have been 
calculated White (8) and Russell (6). 

Line strengths have also been obtained from solar data. The method has 
been fully described Menzel (5) and more recently Wright (9). The 
method obtain the so-called ‘curve growth’ for number lines 
the solar Fraunhofer spectrum, plotting the theoretical strength each 
line against the corresponding measured equivalent width. The equivalent 
width measure the energy removed from the continuous spectrum 
the absorption line. From the ‘curve growth,’ the strength any line 
may obtained from the measured equivalent width. Solar line strengths 
are more numerous than the theoretical values, since all lines that have 
measurable equivalent widths are used. The solar line strengths also include 
intercombination lines, and these are absent from the theoretical values. 


received May 14, 1946. 
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Beginning 1913, and continuing the work over long period years, 
King has estimated the intensities thousands lines the electric 
furnace and arc spectra. description the apparatus used and the 
procedure followed (1) estimating intensities has been given. The eye 
estimates were made remarkably consistent scale, 
effects were indicated and allowed for when present. These data are useful 
for obtaining line strengths. was pointed out Russell (7) 1925 that 
the square King’s estimated intensities are very nearly proportional 
theoretical line strengths. This relation makes possible interpret King’s 
estimated intensities terms line strengths. will indicated later 
that emission line theory yields the same results. 


Calibration the King Intensity Scale 


The elements chosen for the initial work are Fel, Til, and The 
lines are first divided into groups the same excitation potential. For the 
lines each group, logarithmic plot made between the squares King’s 
arc intensities abscissae and the solar line strengths ordinates. The 
solar line strength called Each resulting line segment now differs 
5040 

excitation potential electron volts the lines each group; the 
excitation temperature degrees absolute the arc. The excitation tempera- 
ture the arc determined from the shift required superimpose each line 
segment the segment for the spectral lines lowest excitation potential. 
the shift called then 


the Boltzmann factor, which call the difference 


5040 (2) 


This relation makes possible determine the excitation temperature 
the arc. 

Fig. shows the resulting curves for after superimposing the line 
segments. The points that are far off the line are probably affected 
blending, and the general scatter partly the result departures from 
coupling. The wave-length range for each plot was chosen after studying the 
sensitivity curve the photographic plates used King. The range 
sensitivity the plates used for the wave-length regions chosen small. 
Trial and error methods also show that the best curves are obtained when the 
data are divided this manner. 

Fig. shows the temperature the arc degrees absolute obtained 
from the data the chosen wave-length regions. The temperature seen 
function the excitation potential the upper level involved the 
production each group lines considered, which means that the arc 
not thermodynamic equilibrium. The excitation potential the upper 
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When the temperature the arc known from Fig. obtain values 
line strengths from Fig. given line read off the value correspond- 
ing the square the logarithm King’s intensity, and then add the 
Boltzmann factor obtained from the temperature the arc and the excitation 


2 
Aac Los 


Fic. Calibration the King intensity estimates. 


of Aac 


Fic. Excitation temperature the arc. 


potential the upper level involved the production the line. This 
calibration extremely useful, since can now obtain strengths for many 
hundreds lines the wave-length range 6600. order 
facilitate more rapid use the calibration, double entry table con- 
structed. Tables and I(c) give the value the line strength desig- 
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nated corresponding given arc intensity and given excitation 
potential. Tables and III give the same quantity for Til 
and Till. 

TABLE (a) 
VALUES 
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4800 A5700 


The question arises whether not the approximation (J)? proportional 
line strength sufficiently accurate. The emission spectral line ergs 
per second from emitting volume cm.’ given the following relation: 


the frequency the line, 
the number atoms the ground level for particular stage 
ionization the element question, 
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the statistical weight the ground level, 


the excitation potential the upper level involved the production 


the line, 
the excitation temperature the emitting volume, 


the theoretical line strength. 
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TABLE 


N 


5.8 


5.4 


4.6 


won N 
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LABORATORY STUDY VISIBILITY THROUGH 


Abstract 


Various opinions have expressed regarding the importance such factors 
‘diffusing due small angle scattering determining the visual 
Ifareduction apparent contrast and influence bright- 
ness level are the only important factors, conditions the obscuration point 
should described by, 


where denotes the brightness object which has contrast with its 
background, and represent respectively the cloud brightness and optical 
density the line sight, and denotes the contrast limen Experiments 
test this relation have been performed with ammonium chloride smoke 
results indicate that the equation adequate and 
that factors not taken into consideration play negligible role total 
obscuration under laboratory conditions. experiments were performed 
various brightness levels within the 100 millilambert range. 


The theory the visual range cloud such smoke fog has received 
considerable attention (for summary see reference (7) The presence 
cloud between observer and object influences the brightness level 
which the observer’s eye must operate and results apparent contrast 
between the object and its background that less than the actual value. 
the cloud may operate obscure the object through additional 
factors such the production small angle scattering apparently diffuse 
object boundaries. This ‘diffusing effect’ has been stated (1) play 
important part obscuration but more recent theoretical calculation (8) 
indicates the converse. 

the obscuring effect cloud determined solely the reduction 
apparent contrast and the influence the brightness level, follows that 
the obscuration point, 


where denotes the actual brightness the object (assumed greater than 
that the background, B,), and and denote respectively the cloud 
brightness and optical density the line sight. represents the actual 
contrast between the object and its background, and defined 


(2) 
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The ‘contrast limen’ denotes the maximum contrast that imperceptible 
the function the brightness level and the angular size 
the object, and determined the observer’s visual ability and not 
any characteristic the cloud. the obscuring effect cloud dependent 
other factors Equation (1) must extended the addition appropriate 
terms. ‘The experiments which this report performed 
ascertain whether the addition such extra terms necessary describe 
obscuration data obtained smoke chamber, i.e., whether reduction 
apparent contrast and influence the brightness level are the only important 
factors determining the obscuring effect cloud under laboratory conditions. 
(1) may rewritten 


the object’s subtended angle and the brightness level, and denotes the 
value the limen some selected angle and brightness level. Equation (3) 


where log For brightness levels greater than about mL. 
0 


and subtended angles greater than 0.6°, 0.01 less. Under these 
conditions little error made regarding constant for objects having 
contrast 0.5 more. plot against log the optical density 
smoke required obscure object (with small correction for the 
log term) should this basis yield straight line unit slope and constant 
consistent with acceptable value the limen the underlying 
assumptions are inadequate (i.e., factors other than the reduction apparent 
contrast and the influence brightness level are important) departures from 
these characteristics may expected. Such considerations formed the basis 
for the experiments. 


Apparatus and Procedure 


(a) The Smoke Chamber 


The main features the smoke chamber and its accessory apparatus are 
shown diagrammatically Fig. The dimensions the chamber were 
1.83 1.22 1.22 m*. was painted flat black inside. The observer, 
seated his observations through the glass window W,. The 
brightness the smoke the observer was brought the desired value 
shining light from banks incandescent lamps and into the smoke 
through the glass windows windows were fitted outside metal slides 
the chamber walls that they could cleaned any time during 
experiment without disturbing the smoke. The test object the form 
shown Fig. was enclosed smoke-tight box having glass window 
was illuminated shielded battery-operated lamps and L’. Since 
the test object was illuminated also light scattered from the smoke, provi- 


LANGSTROTH AL.: VISIBILITY THROUGH CLOUDS 


sion was made for measuring the illumination intensity swinging 
photocell into position before section No. when desired. Both the box 
and the window were located centrally the cross-section the chamber. 
The former could moved along the chamber axis means pipe 
which also served conduit for connecting wires. The pipe carried scale 
from which the depth the smoke layer between the observation window 
and the test object could determined. chamber was fitted with 


Fic. chamber and accessory apparatus; lop view, approximately scale. 


Fic. The test object. Sections and were white and section No. was 
uniform grey. The numerals shown were nat included. The brightness contrast between 
sections Nos. and was for the most part 0.65. 


photoelectric parallel-beam device for measuring the extinction 
cient the smoke was introduced through inlet the bottem 
the chamber and removed when desired through outlet with the aid 
centrifugal blower. fan mounted just below the observation window 
served dilute the smoke during production and stir the contents the 
chamber when desired. 


(b) Production Smoke 

The apparatus used for producing ammonium chloride smoke shown 
compressor was controlled valves and and measured the 
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The meter was modification device (2) used for measuring the 
velocity air currents not bounded confining walls. The air was preheated 
exit temperature 300° passing through lagged Pyrex tube J/, 
cm. long, containing internal ohm heating coil. maintain this 
exit temperature the power dissipation was about 430 watts for air flow 
litres per second. Ammonium chloride was placed Pyrex tube 
forming the inner walls the ohm resistance furnace grams was 
decomposed min. with power input 230 The decomposition 
products mixed with preheated air flowed through the lagged Pyrex tube 
cm. long, which had ohm external heating coil sufficient maintain the 


Fic. Ammonium chloride smoke generator; side view, approximately 


temperature While the furnace was heating, preheated air escaped 
through the outlet but when smoke appeared this was closed and the outlet 
into the chamber opened. The mixture decomposition products and hot 
air was rapidly diluted and cooled blast from the fan Fig. blowing 
directly across the outlet The temperature rise within the chamber was 
about during the production cloud high optical density 
3). 


(c) Procedure 


general two types experiment, each consisting long series 
observations, were performed. the first type the cloud brightness 
was maintained constant throughout the series, the test object brightness 
(or the depth the smoke layer before the test object) was set successively 
various selected values, and the depth the smoke layer (or the test object 
brightness) was adjusted each step until the test object was just invisible 
the observer. the second type experiment, was left unaltered 
throughout series measurements except for variations due changes 
the scattered light incident the test object, was set various selected 
values and obscuration was obtained each adjustment judging 
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the obscuration point, observations were made each two criteria, (a) total 
disappearance, and loss definition. For the former, conditions were 
adjusted until the test object was entirely invisible; for the latter they were 
adjusted until section No. the object appeared light patch undefined 
shape. While observers were perfectly adapted the prevailing brightness 
level experiments characterized constant value, there may have 
been slight departures from complete adaptation experiments which 
was varied. 

cycle observations consisting measurements all selected values 
the predetermined quantity was repeated three five times during the life 
smoke, yielding determinations the conditions required for 
obscuration. Since the size and number the smoke particles varied with 
time, determinations the extinction coefficient were made 


Time, min., from arbitrary zero 


Fic. typical extinction coefficient vs. lime curve, 


and the value the time each observation was read from vs. time’ 
curve plotted from the results. typical curve shown Fig. For the 
same reason, frequent adjustments were required the illumination provided 
the bank lamps (G, G’, Fig. give the desired cloud brightness. 
The adjustments were made changing the number lamps each bank 
and, order avoid noticeable colour differences the illumination, only 
minor current variations were used. Adjustment the test object brightness 
was made varying the current through the lamps (Fig. 1). 

measurement cloud brightness was made immediately before each 
obscuration determination with the aid flicker brightness meter (4) 
sighted through the observation window from just behind the 
observer. Similarly reading the test object brightness was made 
swinging the photocell (P, Fig. into position front the test object and 
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reading the response microammeter; the response was converted terms 
the brightness section No. with the aid calibration curve for 
‘response vs. brightness section No. read with flicker brightness meter 
the absence smoke. The brightness contrast between sections Nos. 
and the target was calculated from numerous flicker meter readings 
the absence smoke. 

The experiments provide data the values and obscura- 
tion over wide ranges all these variables. All brightness values were based 
essentially the readings the flicker meter. The relative values bright- 
ness (which are primary importance this work) are error not more 
than absolute values may more greatly error since absolute 
calibration the meter was made against commercial illumination meter 
with the aid magnesium oxide screen. 


Correction for Limen Variations 

Although variations the limen with the angle subtended section 
No. the test object with and with were not excessive these 
experiments, they could not disregarded. The log term (4) 
may regarded correction applied the observed values. 
The size the correction was some cases great 10%, though generally 
was much less (cf. Table I). 

order obtain correction data the observers determined relative limen 
values for various the brightness level range the experiments. 


TABLE 
DETAILED DATA FOR EXPERIMENT 61A* 


t D Q B. 

Uncorr. Corr. 

1.87 0.50 0.92 5.16 0.94 0.95 

1.76 0.76 0.82 44.1 5.23 1.34 1.37 
1.70 0.97 0.76 12.9 5.00 1.65 1.71 
1.60 1.28 0.69 5.20 4.95 2.05 2.18 
1.45 0.65 2.04 3.01 2.20 2.41 
1.46 0.72 0.84 106 6.60 1.05 1.07 
1.40 0.93 0.77 39.8 5.60 1.30 1.35 
1.34 0.69 12.3 1.68 1.79 
1.26 1.45 0.65 4.86 2.05 1.83 2.01 
1.45 0.65 2.00 1.29 1.77 1.99 
1.15 0.89 0.78 1.03 1.06 
0.98 0.76 37.8 4.00 1.09 1.14 
1.06 1.24 0.69 2.61 1.32 1.43 
0.95 1.16 96.2 9.5 1.10 1.16 


was sel selected values and adjusted obscuration each. Corrections the 
values have been made with the aid Fig. time which the observation was 
made was reckoned from arbitrary sero set soon after the production smoke was complete. 
Brightness given millilamberts, distance metres, subtended angle degrees, and extinclion 
coeficient data refer total obscuration. The observer was G.L. 
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sheet glass was placed before the observation window angle 45° 
the line sight. large brightly illuminated white card was placed 
that the observer saw its image superimposed the test object. With 
smoke the chamber, the apparent contrast section No. was reduced 
until was just imperceptible; this was done either increasing the bright- 
ness the veiling glare provided the white card decreasing the 
brightness the test object. The limen was calculated from 
+B’), where denotes the brightness the veiling glare. Determinations 
were made various test-object-to-observer distances using total disappear- 
ance judging the critical point. The average results obtained were 


follows. 
Subtended angle: 0.73° 0.65° 


Relative limen: 1.00 1.03 1.36 
The angle refers that subtended side section No. Data the 
relative limen values various brightness levels were taken from the literature 
(3, 9*, The average values for subtended angles more are 
follows. 

Brightness level, mL.: 100 1.0 
Relative limen: 1.00 1.01 1.04 1.13 1.34 
Combined with the preceding figures for angle variation these yield data 
Equation (4). was taken the limen for brightness level 100 mL. 
and subtended angle 0.99°. The data used making the corrections 


0.3 


\ 


0.2 
0.0 
100 
B., mL. 


Fic. Variation the limen with subtended angle and brightness level used cor- 
observed values. 


are presented Fig. The appropriate log values were added the 
observed values. The absolute values the limen did not enter into the 


considerations. 
Results 


The data for typical experiment are TableI. this experiment, 
all others, the angle subtended the observer side the illuminated 


See paper Blanchard (Phys. 1918) for the brightness units appropriate 
Konig and limen data. 


4 
q 


CANADIAN JOURNAL RESEARCH. VOL. 25, SEC. 


patch smoke (i.e., the observation window) was about surround 
beyond the window was relatively dark. The typical scatter encountered 
these experiments indicated the Table data plotted Fig. 
summary the results six experiments three observers given 
Table 


corrected 


0.03 0.1 1.0 


Fic. typical plot observed data (Expt. 61A, observer illustrating the linearity 
the relation, the scatter observations, and the range covered. The points for each cycle 
measurements are represented separate symbol. 


TABLE 
SUMMARY THE RESULTS SIX EXPERIMENTS* 


Experiment 54A 56A 62A 60A 61A 
range, mL. 16.8 20.6 1.7 -62 1.3 -76 2.0 -111 


A(uD uncorr.) 1.76 1.82 1.79 1.81 1.69 2.04 
(wD uncorr.) 0.82 0.86 0.80 0.76 0.85 0.85 
corr.) 1.93 1.94 1.90 2.04 1.80 
corr.) 0.99 0.95 0.88 0.95 1.00 1.02 
Averages 

1.82 1.99 

0.82 0.97 


The refer total obscuration. They are described the relation log 
whether corrected for limen variation with subtended angle and brightness level 
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Discussion 


The experimental data exhibit moderately large scatter illustrated 
the typical example shown Fig. The scatter largely the result 
momentary fluctuations the visual ability observers (cf. (3) and can 
hardly avoided experiments this kind. the average, however, 
the relation between the various observed quantities well represented 
the data Table II. 

The range covered the data reasonably wide, i.e. brightness level, 
1.3-111 test object brightness, 0.35-34 mL.; extinction coefficient, 
0.72-2.47 optical density, 0.44 2.28; depth smoke, 0.30 
1.45 For total obscuration the object, the results are described 


where has been ‘corrected’ the addition the log term the observed 
values. equation identical form with Equation (4) and the coeffi- 
cient the log term unity within the limits experimental error. 
Furthermore, equating with 0.65 the constant term one 
obtains value 0.0067 for the limen 100 mL. and 0.99°. This 
excellent agreement with the average value for unrestricted binocular vision 
therefore concluded that Equation (1) adequate and that factors such 
‘diffusing effect’ caused small angle scattering play negligible role total 
obscuration under laboratory conditions. 

Observations the apparent loss definition test objects contrasted 
with total disappearance are discussed the following article. 
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THE RECOGNITION OBJECTS NEARLY OBSCURED 
CLOUD! 


Abstract 


the apparent contrast square test object situated behind smoke layer 
decreased laboratory experiments point reached which the object 
appears light patch unrecognizable shape. decreasing the apparent 
contrast little further the object becomes completely invisible. study 
the conditions under which the apparent loss definition occurs was carried 
simultaneously with the study total obscuration described the preceding 

aper. The optical density smoke required produce loss definition was 
ound linearly related the logarithm the ratio object cloud 
brightness, close correspondence with the relation found for total obscuration. 
Over the range our experiments the minimum optical density smoke 
associated with loss definition was from 23% less than that associated 
with total obscuration, the magnitude the difference depending brightness 
conditions. Since the phenomenon loss definition observable the 
presence veiling glare without smoke, probably physiological origin. 


Experiments the total obscuration square test object ammonium 
chloride smoke determined extinction coefficient have been described the 
preceding paper. the experiments the layer depth the brightness 
the smoke the illumination incident the object was adjusted until the 
object became just invisible. Information was obtained the brightness 
conditions and optical density smoke required for complete obscuration 
over fairly wide range the variables involved. the same experiments 
analogous data were obtained the conditions associated with apparent loss 
definition. The adjustments were made until the shape the object 
became just unrecognizable although its presence was still indicated light 
patch undefined shape. The apparatus and procedure are described the 
preceding article. the purpose this note place the data record. 


Results 


Data for typical experiment are presented Table and Fig. Table 
and Fig. the preceding paper contain the corresponding data for complete 
obscuration the saine experiment. Table II, corresponding Table 
the preceding paper, contains summary the results seven experiments 
three observers. The symbols and have their previous signi- 
ficance. denotes the minimum depth smoke layer required produce 
loss definition. 
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TABLE 


DETAILED DATA FOR EXPERIMENT 61A* 


1.91 0.40 0.95 89.1 5.40 0.76 
1.80 0.60 0.87 43.6 5.70 1.08 
1.71 0.82 0.80 13.2 5.28 1.40 
1.63 0.72 5.00 5.40 1.81 
1.37 0.67 2.06 5.55 2.14 
1.49 0.60 0.87 101 6.12 0.89 
1.41 0.80 0.91 40.8 5.68 1.13 
1.35 1.10 0.73 12.6 5.28 1.48 
1.27 1.36 0.67 5.28 5.49 1.73 
1.24 1.45 0.65 2.00 2.98 1.80 
1.18 0.64 0.85 92.9 4.27 0.76 
1.13 0.81 0.80 37.2 2.43 0.92 
1.09 1.06 0.73 13.8 2.86 1.16 
0.96 0.87 0.79 94.5 7.60 0.83 


was set selected and adjusted the point loss definition each. The 
(min.) which the observations were made were reckoned from arbitrary zero set soon 
after production smoke was complete. Brightness given distance metres, 
subtended angle (a) degrees, and extinction The observer was 
The total obscuration data for this experiment are given Table the preceding paper. 


0.1 1.0 
(B./B.) 


Fic. typical plot observed data (Expt. 61A, observer G.L.). The points for each 
cycle measurements are represented separate symbol. 


Discussion 


Fig. provides The scatter plotted points comparable 
that obtained total obscuration determinations (compare Fig. with 
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TABLE 


A SUMMARY OF THE RESULTS OF SEVEN EXPERIMENTS* 


Experiment 55A 62A 60A 61A 
range, mL. 16.8 17.5 20.6 1.7-62 1.3-76 12.3 


0.50-24 0.98-34 0.77-15 0.50-7.0 0.35-6.3 1,3-9.5 


A’ 1.55 1.55 1.62 1.56 1.37 1.77 1.71 
E’ 0.77 0.83 0.80 0.70 0.71 0.77 0.78 
Averages A’ = 1,59 E’ =0.77 


The data plotted for each experiment were found well described the relation 
B.); see Fig. The data total obscuration for the above experi- 
ments may found Table the preceding paper. 


Fig. the preceding paper). Using the average results given Table II, 
the minimum optical density required produce loss definition 


described 


This equation simply statement average experimental results; 
correction corresponding the correction for limen variation applied 
the preceding paper, has been made. The equation therefore 
includes any variations originating changes subtended angle (about 
0.6° 1°) and brightness level (about 110 


comparison between the optical density required for loss definition 
and that required for complete obscuration readily made. The latter 
described Equation (1) with 1.59, and 0.77 replaced respectively 
1.82, and 0.82 (Table the preceding paper). The values are 
uncorrected for limen variation. Subtraction the two equations yields, 


where denotes the difference between the optical density required 
produce complete obscuration and that required for loss definition under 
the particular conditions our experiments. the lowest and highest 
ratios used (0.044 and 3.16) the values were respectively and 
12% less than the values. Any effect produced variation subtended 
angle the test object contained implicitly Equation (2). 


described the preceding paper, corrections for limen variations with 
subtended angle and brightness level were applied the total obscuration 
data. corresponding corrections have been attempted for the loss 
definition data because the lack appropriate information. determin- 
ing limen values with the aid veiling glare the absence smoke, how- 
ever, was noted that loss definition preceded total disappearance the 
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apparent contrast the object was gradually decreased. The corners the 
square test object became indistinct that might easily have been confused 
with disk although its presence was not doubt. Since the phenomenon 
loss definition was observed the absence smoke, considered 
physiological origin; probably not connected with any objective 
factor such diffusing effect caused the smoke. 
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NOTE THE DETERMINATION THE EXTINCTION 
COEFFICIENT 


Abstract 


Difficulties were encountered smoke chamber experiments obtaining 
reliable extinction coefficient data methods based brightness measure- 
ments bright and ‘black’ surfaces seen through the smoke. These difficulties 
led investigation the reliability commonly used method employing 
parallel light beam and photoelectric cell. The results indicate that the 
method capable yielding extinction coefficient data with accuracy 
under the conditions our experiments. 


surface brightness seen through smoke layer depth and 
extinction coefficient has apparent brightness given 


where denotes the smoke brightness the line sight. The extinction 
coefficient described 


and can determined the other quantities are measured. attempt 
use this method smoke chamber experiments with smokes having coefficients 
the range 0.7 2.5 proved unsuccessful. The cm.? surface 
was located about cm. behind the glass window small smoke-tight box 
that could moved along the chamber axis alter the d-value, and was 
illuminated shielded lamps within the box. The smoke chamber 
described elsewhere (1). Brightness measurements were made with flicker 
brightness meter (2). Determinations made with the meter sighted 
just beyond the boundary the surface led which depended 
and were obviously incorrect. When the determinations were made 
with the meter sighted the opening cm. diameter tube located 
near the centre the surface and serving ‘black’ body, the resulting 
were acceptable this respect. They were however 30% 
lower than those expected from the slope curves for ‘log vs. layer depth 
obscuration’ obtained with constant veiling glare produced the 
illuminated white card 45° glass sheet front the 
observer. These difficulties were associated with inability make 
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appropriate determinations cloud brightness, and led investigation 
method employing parallel beam light and photocell. The investi- 
gation forms the basis this report. 


The Nature the Experiments 


(a) The photoelectric parallel-beam method.—This method, which not 
original with us, has been used several laboratories, but not know 
the name the originator. The experimental arrangement used 
shown diagrammatically Fig. collimating and condensing lens 


Fic. the experimental arrangement for the parallel-beam method extinction 
coefficient measurement. 


and L’, focal length cm.) were placed the opposite sides the smoke 
chamber, distance 1.30 apart. cp. automobile lamp placed 
the principal focus the lens and operated constant voltage served 
the light source screen containing small opening was placed the 
principal focus the lens The purpose the screen was prevent 
much the scattered light possible from reaching the vacuum photocell 
while permitting the directly transmitted light forming the image the 
filament so. The relation between the amplified photoelectric 
current and the light energy incident the cell was determined over the 
required range with the aid large Nicol prisms inserted between and L’. 
Denoting the light energy incident the cell the presence and the 
absence smoke and respectively, the extinction coefficient the 
smoke was calculated from 

(b) Dependence the width the collimated light order 
for the method prove satisfactory the contribution directly transmitted 
light the energy incident the photocell must predominate. Compared 
the contribution light that has undergone scattering within the beam, 
has been scattered back from regions outside, must negligibly small. 
The directly transmitted light forms image the lamp filament the 
plane the screen but the scattered light incapable doing so. The 
contribution scattered light the total incident the photocell expected 
depend the beam width (Fig. 1), approaching zero for infinitely 
narrow beam. observed dependence beam width might therefore 
taken indication that scattered light cannot neglected. Experi- 
ments test this point were performed with the aid diaphragms various 
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Sketch the surface distribution light incident photographic plate placed 
the plane the screen Theimage formed directly transmitted light shaded, 
and the incidence scattered indicated the dotted area. The dashed circle indicates 
the position the opening the screen. The plate was photometered along the band AB. 


100 


Relative intensity 


Distance along plate 


the main image the intensity low, and measurements elsew have been 
omitted avoid confusion. The unplotted points fall closely along those plotied. The dashed 
line suggests upper limit for the intensity scattered light the area occupied the 
filament image. 
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sizes inserted the collimating lens. Typical results are given Table 
Section 

(c) Estimation upper limit for the contribution scattered 
order obtain estimate the relative energies associated with directly 
transmitted and scattered light admitted the photocell, the screen (Fig. 
was replaced photographic plate. short exposure was made the 
early stages smoke using beam width Calibration marks 
were placed the plate with the aid adjustable high speed sector and 
light source appropriate distribution. These permitted 
photographic densities the plate converted terms the relative 
light energies incident. The plate (Eastman 40) was somewhat under- 
developed produce low contrast, care being taken minimize the Eberhardt 
effect. Microphotometric determinations the photographic density were 
made intervals 0.28 mm. along the band indicated Fig. The 
determined distribution light energy along the band Fig. shown 

(d) The dependence the angle divergence non-parallel 
permissible, the use somewhat divergent beam light has certain advan- 
tages since the adjustments are less critical and the apparatus consequently 
more rugged. Series measurements were made with beams various 
divergence. The divergence was controlled adjusting the distance between 
the light source and the collimating lens. Typical results are given Table 


Section 
Results 


TABLE 


A COMPARISON OF THE M-VALUES OBTAINED WITH PARALLEL BEAMS OF VARIOUS 
WIDTHS IN TYPICAL EXPERIMENTS* 


for beam widths 


3.94 3.30 2.06 1.31 0.88 
8.5 1.82 1.80 1.90 1.95 1.84 
16- 1.87 1.85 1.90 1.93- 1.88 
1.74 1.72 1.79 1.75 1.68 
71- 1.67 1.64 1.69 1.64 1.65 
106-110 1.48 1.45 1.49 1.50 1.47 
49A 12- 2.39 2.38 2.38 2.33 
65- 1.77 1.76 1.76 1.76 
48A 3.5-6 2.81 2.88 2.95 2.88 
40- 42.5 2.80 2.86 2.87 2.85 
111-112.5 2.37 2.42 2.45 2.43 


Beam width given cm., and denotes the time interval which the observations 
were made reckoned from arbitrary zero set soon after complete generation ammonium 


chloride smoke. 


Experi- 
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TABLE 


COMPARISON OBTAINED WITH BEAMS VARIOUS ANGLES DIVERGENCE* 


for divergence angles (degrees) 


0.0 2.6 4.6 7.6 10.0 15.6 20.0 30.0 
13-19 1.85 1.85 1.84 1.71 1.67 1.63 1.54 
30-36 1.63 1.64 1.62 1.55 1.53 1.50 1.43 
53-59 1.47 1.47 1.43 1.34 1.31 1.22 1.19 


The angle referred the plane angle the apex the cone back- 
ward the boundaries the ray bundle between the collimating and condensing lenses. The time 
interval has the same significance Table For values 4.6° less the collimating lens 
down diameter for the others the diameter the diaphragm was 


Discussion 


The particulate characteristics smoke change with time. the 
experiments each group measurements compared was made rapidly 
possible minimize this factor. Examination the decrease extinction 
coefficient with time for given beam width (Table indicates that changes 
smoke characteristics during the time interval required for group 
comparative determinations caused decrease not more than 0.01 0.02 
the extinction coefficient value. 

The results Table show systematic variation the determined 
extinction coefficient value with beam width. The average deviation 
individual from the means for the groups 1.2%. This evidence 
supports the view that the contribution scattered light the energy incident 
the photocell negligible the parallel-beam method. 


The contour Fig. indicates the surface distribution energy incident 
the screen before the photocell. probable upper limit for the scattered 
energy the region the filament image indicated the dashed line 
the figure. comparison the energies scattered and directly transmitted 
light admitted the photocell was made integration under the respective 
contours over the area the screen opening. The result indicates that the 
former was less than the latter. This value upper limit. The 
occurrence halation diffusion light traversing the emulsion are 
factors that would influence the extrapolation the scattered light contour 
toward spuriously high values, and the use colour blind emulsion probably 
favoured the scattered the expense the transmitted light. The Eberhardt 
effect would operate the opposite direction, but evidence its presence 
energy introduces error extinction coefficient 2.00 and 
error coefficient 1.00. The result indicates that the parallel- 
beam method capable yielding accurate values for the extinction co- 


efficient. 
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The results quoted Table show that not important use strictly 
parallel beam under the conditions these experiments. Divergences 
much appear permissible. 

Extinction coefficients measured the parallel-beam method were found 
about 20% higher than those determined with the aid brightness 
meter, illuminated surface, and ‘black’ body described the introduc- 
tion. mentioned there, this finding accord with the results obscura- 
tion data. 

References 
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Fic. Sketch the surface distribution light incident photographic plate placed 
the plane the screen (Fig.1). Theimage formed directly transmitted light shaded, 
and the incidence scattered light indicated the dotted area. The dashed circle indicates 
the position the opening the screen. The plate was photometered along the band AB. 


100 


Relative intensity 


0 
Distance along plate 


Fic. contour obtained photometering two photographic plates. Outside 
the main image the intensity low, and measurements elsewhere been 
omitted avoid confusion. The points fall closely along those plotted. The dashed 
line suggests upper limit for the intensity scattered light the area occupied the 
image. 
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sizes inserted the collimating lens. Typical results are given Table 
Section 

(c) Estimation upper limit for the contribution scattered 
order obtain estimate the relative energies associated with directly 
transmitted and scattered light admitted the photocell, the screen (Fig. 
was replaced photographic plate. short exposure was made the 
early stages smoke using beam width 3.9 cm. Calibration marks 
were placed the plate with the aid adjustable high speed sector and 
light source appropriate distribution. These permitted 
photographic densities the plate converted terms the relative 
light energies incident. The plate (Eastman 40) was somewhat under- 
developed produce low contrast, care being taken minimize the Eberhardt 
effect. Microphotometric determinations the photographic density were 
made intervals 0.28 mm. along the band indicated Fig. The 
determined distribution light energy along the band Fig. shown 

(d) The dependence the angle divergence non-parallel 
permissible, the use somewhat divergent beam light has certain advan- 
tages since the adjustments are less critical and the apparatus consequently 
more rugged. Series measurements were made with beams various 
divergence. The divergence was controlled adjusting the distance between 
the light source and the collimating lens. Typical results are given Table 


Section 
Results 


TABLE 


A COMPARISON OF THE J-VALUES OBTAINED WITH PARALLEL BEAMS OF VARIOUS 
WIDTHS IN TYPICAL EXPERIMENTS* 


for beam widths 


Experi- 
ment 3.94 3.30 2.06 1.31 0.88 
47A 1.82 1.80 1.90 1.95 1.84 
16- 1.87 1.85 1.90 1.93- 1.88 
52- 1.74 1.72 1.79 1.75 1.68 
71- 1.67 1.64 1.69 1.64 1.65 
106-110 1.48 1.45 1.49 1.50 1.47 
49A 12- 2.39 2.38 2.38 2.33 
65- 1.77 1.76 1.76 1.76 
48A 3.5-6 2.81 2.88 2.95 2.88 
42.5 2.80 2.86 2.87 2.85 
111-112.5 2.37 2.42 2.45 2.43 


Beam width given cm., and denotes the time interval which the observations 
were made reckoned from arbitrary zero set soon after complete generation ammonium 
chloride smoke. 
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TABLE 


COMPARISON OBTAINED WITH BEAMS VARIOUS ANGLES DIVERGENCE* 


for divergence angles (degrees) 
0.0 2.6 4.6 7.6 10.0 15.6 20.0 


30-36 1.63 1.64 1.62 1.53 1.44 1.43 
53-59 1.47 1.47 1.34 1.31 1.19 


The angle referred the plane angle the apex cone back- 

ward the boundaries the ray bundle between the collimating and condensing lenses. The time 

interval has the same significance For values 4.6° less the collimating lens 

down diameter 1.31 cm.; for the others the diameter the diaphragm was 
cm. 


Discussion 


The particulate characteristics smoke change with time. the 
experiments each group measurements compared was made rapidly 
possible minimize this factor. Examination the decrease extinction 
coefficient with time for given beam width (Table indicates that changes 
smoke characteristics during the time interval required for group 
comparative determinations caused decrease not more than 0.01 0.02 
the extinction coefficient value. 

The results Table show systematic variation the determined 
extinction coefficient value with beam width. The average deviation 
individual from the means for the groups 1.2%. This evidence 
supports the view that the contribution scattered light the energy incident 
the photocell negligible the parallel-beam method. 


The contour Fig. indicates the surface distribution energy incident 
the screen before the photocell. probable upper limit for the scattered 
energy the region the filament image indicated the dashed line 
the figure. comparison the energies scattered and directly transmitted 
light admitted the photocell was made integration under the respective 
contours over the area the screen opening. The result indicates that the 
former was less than the latter. This value upper limit. The 
occurrence halation and diffusion light traversing the emulsion are 
factors that would influence the extrapolation the scattered light contour 
toward spuriously high values, and the use colour blind emulsion probably 
favoured the scattered the expense the transmitted light. The Eberhardt 
effect would operate the opposite direction, but evidence its presence 
energy introduces error extinction coefficient 2.00 and 
error coefficient 1.00. result indicates that the parallel- 
beam method capable yielding accurate values for the extinction co- 
efficient. 
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The results quoted Table show that not important use strictly 
parallel beam under the conditions these experiments. Divergences 
much appear permissible. 

Extinction coefficients measured the parallel-beam method were found 
about 20% higher than those determined with the aid brightness 
meter, illuminated surface, and ‘black’ body described the introduc- 
tion. mentioned there, this finding accord with the results obscura- 


tion data. 
References 
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NOTE PHASE CORRECTION ELECTRICAL 
DELAY 


Abstract 


The standard method correcting the phase characteristic low-pass net- 
works introducing mutual inductance between the coils extended. 
means mutual inductance between the coils and capacity across them, 
possible eliminate the third and fifth order terms the variation the 
shift with frequency. The relation required between the circuit constants 
this result given, and found very simple. 


There are number systems that may used provide small delay 
for video signal. These may divided roughly into two categories; 
those with distributed constants, and those with lumped constants. 
the first category are concentric lines with coiled inner conductors and the 
‘condensed cables’ suggested Kallmann (2). These are simple devices, 
fairly readily constructed, and they provide reasonably good solution the 
problem. Their performance not perfect—they suffer phase defects 
frequencies for which the wave-length along the structure comparable with 
its diameter. the second category are wave filters and standard artificial 
lines. with this type that this report concerned. Holcomb 
(1) has described elaborate bridged-T delay network with excellent phase 
characteristics, which, however, much too complicated for high frequency 
ranges. 

well known (3, 245) that the phase characteristic simple low- 
pass filter can improved introducing mutual inductance between adjacent 
coils shown Fig. the phase shift single section and the 
Maclaurin series, obtain 


are the Maclaurin coefficients and are functions the network 
parameters. may choose the mutual inductance that vanishes, 
removing the third order term The condition for this is* 


(2) 


Some attention must given the sign Equation (2) implies that 
the coils are series-aiding. 
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recommends which over-corrected the sense taken here. should noted that 
chain m-type sections would produce network with the inductances divided halves the 
manner suggested the end the paper. 
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will shown how introduce another circuit parameter obtain yet 
another order phase correction, the word ‘correction’ here having the 
significance used above, i.e., the removal successive non-linear terms 
Equation (1). The circuit that Fig. the added element being con- 
denser across the inductance and may chosen that the 
and for this circuit both vanish. The required conditions reduce very 


simple form, namely 6.439 (3) 


Here again found that the coils are series-aiding. 


Fic. low-pass filter delay network with mutual inductance for phase correction. 


Fic. Low-pass delay network using both mutual inductance and capacity, Ci, for 
phase correction. 


This result obtained laborious but straightforward algebra. Only 
outline the calculation will given. network equations are: 


put 


where and are arbitrary constants and the propagation constant, 


then the consistency Equations (5) and (6) requires that 


Equation (7) determines the propagation constant 
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there are two solutions for cosh there are two modes transmission. 
One these rapidly attenuated all frequencies and does not concern us. 
The other mode has low-pass filter type behaviour that susceptible 
the improvement obtained choosing the circuit constants properly. 

the pass band, purely imaginary and equal The solution 
Equation (7) then gives cos This solution may expanded 
Maclaurin series 


The coefficients are all functions the network parameters. 


may regarded additional arbitrary parameter. 

this point could apply Taylor’s theorem Equation (7) (8) 
obtain expansion the form (1), and then equate and zero get 
our desired conditions. However, much easier proceed follows. 


have 


The required result obtained identifying the coefficients Equations 
(8) and (9) and including the sixth degree terms. That is, put 


Eliminating between these equations gives two equations for the circuit 
parameters, These can reduced 


6.439 


28.439 
that is, Equations (3) and (4). 
Turning other important constants the circuit, may readily find 


that the cut-off frequency 


The question the characteristic impedance the proper terminating 
impedance more difficult. the mutually coupled coils replaced 
their equivalents, will seen that the ladder becomes not chain four- 
terminal networks, but six-terminal networks. this fact that 
responsible for the two modes propagation. Such network requires two 
characteristic impedances. 
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practical case the line will probably broken off one the junction 
points they appear Fig. This procedure destroys the six-terminal 
character the last section, but there does not seem any good reason 
for attempting preserve it. The terminating admittance which will 
remove the reflected propagated wave for such arrangement can cal- 
culated. The result appears complicated—it has not been carried 
through any case. will undoubtedly mainly resistive, value close 
Experience with networks with cut-off frequency mega- 
cycles per second indicates that the results are uncritical the termination. 


has been assumed that there mutual inductance between any pairs 
coils except adjacent ones. the coils are mounted coaxially, this will 
not true, i.e. there will interaction between non-adjacent coils. The 
effect this has not been considered, but will probably small. any 
concern were felt about this point, possible divide each coil half and 
then arrange the parts that half-coil one section interacts only with 
the adjacent half-coil the next section. 
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